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Topological aspects of surface states in semiconductors are studied by an adiabatic deforma- 
tion which connects a realistic system and a decoupled covalent-bond model. Two topological 
invariants are focused. One is a quantized Berry phase, and the other is a number of the edge 
states. A winding number as another topological invariant is also considered. The surface states 
of Si and Ge at (111), (HO), and (100) surfaces are classified by the topological invariants. Sur- 
face states of the GaAs as heterosemiconductors are also discussed. 

KEYWORDS: semiconductor surface, Berry phase, edge state, dangling bond, Ge, Si, GaAs 



1. Introduction 

For a few decades the topological feature which does 
not depend on the detail of systems has been attracting 
attention in the condensed matter physics. One of most 
remarkable examples is a quantum Hall effect (QHE), 1 
where a topological invariant known as the Chern num- 
ber 2 has been observed experimentally. The Chern num- 
bers are denned for a bulk state without boundaries. In 
the QHE, the edge states which are characteristic for a 
system with boundary also have a significant topological 
importance. 3 

Edge states or surface states, which localize at a 
boundary of a system, have significant effects in phys- 
ical phenomena; e.g., the zero bias conductance peak 
structure in anisotropic superconductivity, 4, 5 single-layer 
graphite, 6 and the boundary local moment at an edge of 
hexagonally bonded honeycomb sheets consisting of B, 
N, and C atoms. 7 Especially, in a graphene, it is now 
widely understood that there exist important edge states 
which are related to the chiral symmetry of the system 
and the Dirac cone spectrum. In metal-semiconductor in- 
terfaces, an interface state, known as the metal induced 
gap state, has an important role in understanding Fermi- 
level pinning. 8 

Moreover, in realistic materials, there exists a surface 
reconstruction or a surface relaxation, which is sensitive 
to the detail of systems. However, we shall discuss that 
topological invariants which is independent of the detail 
of models is useful for understanding the surface states of 
3 dimensional semiconductors. A candidate of topologi- 
cal invariants for semiconductors is the quantized Berry 
phase which has been proposed as a local order param- 
eter for clarification of the phase of the gapped quan- 
tum liquids. It has been demonstrated that quantization 
of this local order parameter has an advantage over the 
usual order parameter in a frustrated spin system, 9 and 
a strongly correlated system. 10 

In this paper, the quantized Berry phase is evaluated 
for surface states in typical semiconductors such as Si, 
Ge, and GaAs. We note that the Berry (or Zak) phase 
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also appears in the King-Smith- Vanderbilt formula in the 
theory of macroscopic polarization. 11 ' 12 In our case the 
Berry phase is defined for a 1-dimensional system with 
2-dimensional wave vectors to characterize the surface, 
and quantized due to a anti-unitary symmetry. In addi- 
tion, we focused on the number of the edge states and a 
winding number as the other topological invariants, fol- 
lowing Ref.13. For simplification, we limit ourselves to 
surface states in ideal surfaces understood as Shockley 
states. 14 Shockley states at a metal-semiconductor con- 
tact are closely-linked to dangling bonds. 15 Although we 
does not consider a surface reconstruction, the study on 
surface states at ideal surfaces sheds a light on real mate- 
rials in the sense that the number of dangling bonds in a 
respective surface acts as a trigger for surface-dependent 
physics; for example, (7 x 7) structures on (111) surface 
of Si 16 and Ge, 17 (16 x 2) structures on (110) surface of 
Si 18 and Ge, 19 and c(4 x 2) structures on (100) surface 
of Si 20 and Ge. 21 

2. Model and Results 

To consider (111), (110), and (100) surfaces of Si, Ge, 
and GaAs as typical IV and III-V semiconductors, we 
use tight binding models determined by the Slater and 
Roster's energy integrals. 22 We also consider a spin-orbit 
coupling that is important for Ge. We use the tight bind- 
ing and spin-orbit coupling parameters of Grosso and 
Piermarocchi 23 for Ge, those of Klimeck et. al. 24 for Si, 
and those of Boykin et. al. 25 for GaAs. Note that a sp 3 
model is considered for Ge, and a sp 3 s* model for Si, and 
GaAs. 

To define a surface-dependent Berry phase, let us start 
from the definition of a bulk Hamiltonian. A position of a 
center of a unit cell is given by = Y] _-, i n a n , i n G Z, 
where a%, ai, and are primitive vectors. Primitive re- 
ciprocal lattice vectors b n are defined by a n -b rn = 27r5 nm , 
where S nm is the Kronecker delta. A wavenumber vector 



k h 



The bulk Hamiltonian in 



k is given by k = J2' r 
three dimensional periodic boundary condition (PBC) 
Hpbc is given by 

HpbC = C l/3 [ H PBc(k)]p jl3 , C h ,0>- (1) 

k /3,f3' 

Here, the label [3 = (a, a) distinguishes spin-dependent 
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orbitals in a unit cell. We suppose two atoms in a unit 
cell and denote orbitals as a = 1, 1, 2, 2, . . . , N, N, where 
N is a number of orbitals in a atom and a on a atom 
means an orbital opposite to an orbital a on the other 
atom (See Fig. 1). Then, Hp B c(k) is a AN x AN matrix. 




Fig. 1. Model of a diamond structure. Spheres stand for atoms, 
and lines stand for sp 3 hybridized orbitals. A shade stands for 
(111) surface, and a n stands for basis lattice vectors relative to 
(111) surface. 



Next, we suppose that cti and a<i span a plane parallel 
to a surface of a semiconductor as a shaded plane shown 
in Fig. 1. 2-dimcnsional wavenumber vector is given by 
k|| = J2n=i knb n - Then, the bulk Hamiltonian is written 

as 



WpbC = S£4,|,n [#PBc(fc||)L n , C fe|| 



(2) 



where the label n is defined as n = (23, (3) and HpBc(k\\ ) 
is a ANL x ANL matrix, and i 3 e [l,L]. L is a num- 
ber of layers. To consider the surface, we introduce the 
open boundary condition (OBC) for the one dimensional 
system in the a 3 direction. The matrix i?oBc(fe||) is 
given by truncating 7JpBc(fe||)- For example, a natu- 
ral way to truncation is to prohibit all the matrix ele- 
ments across L. We define a boundary-hopping matrix 
iJ B (fcj|) = i/pBc(fc||) - #OBc(&||) as 



[#B(fc||)], 



= / [tfPBc(fe||)] n , n , 





[*3,*3] 3 L 

others 



(3) 



-ff B (fc||) is ANL x ANL matrix. This matrix represents 
all hopping elements across the boundary of the system. 

To define the Berry phase, we introduce a twist angle 9 
in a hopping term across the boundary as c£ n Ck 9 ,n' — * 

e'^cj, n Cfc |hn '. In detail, the twist angle 9 is introduced 

in a selected element nn' of J?b (fell ) and we denote it as 
Hq(9 7 feii). Then, the Berry phase can be defined as 

7n„'(*||)= / (gs\d e \gs)de, (A) 
Jo 

where \gs(9)) is the half-filled ground state of 
-flpBc(#,fc||) = -ffoBc(fe||) + H B (0,k»). Especially, we 



denote main elements of the Berry phase as j aa — 
l(Lacr)(iaa)- The Berry phase is quantized because the 
one-dimensional system at fe|| defined by £/pBc(fe[|) has 
a inversion symmetry at the boundary: (23, a) <-» (L + l — 
i3,a). When we denote the inversion as U and the com- 
plex conjugate as K, the anti-unitary operator 6 = K U 
is commutable with the Hamiltonian. Then, the Berry 
phase turns out to be quantized in the same way as in 
Ref. 9 

Figure 2 (a) shows the band structure of 7Ypbc of Ge 
and reproduce the indirect gap. 26 After introducing the 
(111) surface, the band diagram of Hqbc has edge states 
in the band gap, which are doubly-degenerated per spin, 
as shown in Fig. 2 (b). The Berry phase of Ge (111) 
surface is actually quantized and does not depend on fe|| . 
It turns out to be 711 = 71J. = ir on the strongest hopping 
in Hb and zero on the others. We note that the number 
of candidates of 7„ n / is 48 for the parameters we used. 



(a) 




Fig. 2. (a) A band structure of bulk Ge. (b) Energy spectrum of 
Ge (111) surface. The white area denotes band gap of bulk Ge 
and lines do edge states in Ge (111). The notations of points in 
relative surface Brillouin zone are according to Ivanov et. al. 27 
The number of the degeneracy of the edge states is four including 
spin degrees of freedom. 



The quantized Berry phase and the number of edge 
states n e are topological invariant in the following adia- 
batic transformation. We first modify the spin-orbit cou- 
pling term adiabatically to zero preserving the energy 
gap open, that is, without changing the topological in- 
variants. At this stage the Hamiltonian is spin decoupled 
as H = 7i| © Hi . Now we ignore the spin index below for 
simplification. Next, we can modify the hopping terms 
between different orbitals and obtain a Hamiltonian as 
the sum of two band Hamiltonians H = 7Y (1) 0- • -eW (Ar) , 
where TL^ involves hopping terms between a and a 
only. For this process, it is also possible to hold the 
gap open. After the Fourier transformation of the 03 
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direction, H PBC can be written as a 2 x 2 matrix 
Hp^ c (k = (kpk 3 )). We can expand H^ c by the Pauli 
matrices <x by 

H^ c (k) = fl^(fc) ■ a + R { a \k)a , (5) 

where Co is the 2x2 identity matrix. This equation 
indicates a one to one mapping between H^ c and a four 
dimensional vector (Rq,R). Therefore, for fixed fen, a 
loop C a in R 4 space is given as £3 varies over S 1 = [0, 1). 
We note that the four dimensional vectors at k 3 = 
and kz = 1 are the same. Mathematical representation is 
given by 

C : k 3 e S* 1 -» (Rq,R) G M 4 . (6) 

Thirdly, we can modify i? (fc||,fc3) to zero preserving 
the gap open. For the arbitrary surface of Si or Ge, 
Rff\k\\, £3) is zero. The condition that C a is restricted 
in a two-dimensional plane is called as a chiral symme- 
try. 13 The loop is given by 

C-.heS 1 ^(R x ,R y )eR 2 . (7) 

For a 2-dimensional loop, we can define a winding num- 
ber as in Ref. 13. This is the third topological invariant 
we considered in the paper. The winding number W a {k\\ ) 
is defined as the total number of times that the loop trav- 
els counterclockwise around the origin in the xy plane, 
and determines whether H^ c (k\\) has edge states or 
not. Figure 3 shows C a (k\\) of Ge (111) surface. In each 
plane vertical to the z axis of Fig. 3, four loops with fixed 
A; || are displayed and one of these loops encloses the ori- 
gin, R x = R y = at each k\\, i.e., the winding number 
of a link is Wi(fc||) = 1 for arbitrary fe||, which clarifies 
existence of zero energy edge states. The other loops do 
not enclose for each fe||, i.e., W Q=2! 3.4(fe||) = 0, which 
is consistent with absence of zero energy edge states for 
each fen . Note that zero energy edge states are doubly de- 
generated (left and right). This result from the winding 
number corresponds to the number of edge states n e = 2 
per spin. 

Finally, we modify hopping parameters except for hop- 
ping between face-to-face orbitals which locate on near- 
est neighbor atoms to zero adiabatically. We call this 
decoupled Hamiltonian as a covalent-bond model TLcb- 
Edge states of Hcb with the OBC corresponds to dan- 
gling bonds, where dangling bonds are defined as the 
covalent bonds which are cut by the surface. Of course, 
topological invariants 7 and n e are unchanged between 
the adiabatic deformation. Then, j nn / — it implies that a 
covalent-bond bond exists on the link nn' and the integer 
n e corresponds to the number of dangling bonds. 

3. Summary of the other results 

Here, we summarize results of the other surfaces of 
Si and Ge studied in the same manner as Ge (111). As 
shown in Figs. 4 and 5 in addition to in Fig. 2 (b), each 
energy band diagram of edge states of Ge (110), Ge (100), 
Si (111), Si (110), and Si(100) shows different energy 
dependence. These edge states are all doubly degenerated 
due to two sides of surface. In Ge (110), Si (110), Ge 
(100), and Si (100), value of energy split between two 




Fig. 3. Two-dimensional loops £ a (fc||) of the Hamiltonian with 
the chiral symmetry for the Ge (111) surface. The vertical axis 
shows ku and the other axes show xy elements of the loops. 
The vertical axis equals to a part of horizontal axis in Fig. 2 
(b). Vertical line between two solid circles denotes the origin for 
each fcy . Each loops with fixed fey locates on plane parallel to xy 
plane. Each plane has four loops related to four Hamiltonians. 
On each plane, one loop encloses the origin and the others do 
not. 



(a) Ge (110) 




r x 1 mx r 



(b) Ge (100) 




Fig. 4. Energy spectrums of Ge that have ideal surface (110) and 
(100). The notations of points in relative surface Brillouin zone 
are according to Ivanov et. al. 27 The number in the figures is the 
degeneracy of the edge states are eight including spin degrees of 
freedom. At (100) surface, edge states disappear for special values 
of fc|| . 



kind of edge states is roughly estimated as bonding states 
and anti-bonding states of two dangling bonds on the 
surface. Especially, the edge states merged into the bulk 
bands for (100) cases and we cannot define the number 
of edge states n e . 

The Berry phase 7 turns out to be quantized by the 
anti-unitary symmetry and succeed in clarification of the 
three surfaces for Si, Ge, and the decoupled covalent- 
bond model after the adiabatic transformation while n e 
fails for the (100) surface, as summarized in Table I. Since 
it is topological quantity, the result does not depend on 
the detail of the tight binding parameters unless the gap 
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(a) Si (111) 

3 i r 




r M K 

(b) Si (110) 




r X' mx 
(c) Si (100) 




Fig. 5. Energy spectrums of Si that have ideal surface (111), 
(110), and (100). At (111) surface and (110) surface, four and 
eight surface states appear for each k» including spin degrees of 
freedom, relatively. At (100) surface, eight surface states appear 
including spin degrees of freedom for special values of kn. 



remains open. Then, Ge and Si with the same surface are 
equivalent from a viewpoint of the Berry phases. 

After the adiabatic deformation described above, the 
winding number W for the chiral symmetric Hamiltonian 
classifies the different surfaces. The results are also sum- 
marized in Table I. The sign of W determines whether 
the loop C is clockwise or anti-clockwise, and becomes 
important to classify (110) and (001) surfaces. The lat- 
ter is identical to (100) surface. We note that n e can 
be defined for the chiral symmetric Hamiltonian of Gc 
(100) and Si (100), because edge states emerges from the 
bulk bands into the gap in the adiabatic transformation. 
However, we cannot allow such adiabatic transformation 
if we consider topological invariant. 

The decoupled covalent-bond model clarifies the mean- 
ing of topological invariants because edge states are iden- 
tified as dangling bonds. Then, n e of each materials refers 
the number of dangling bonds, and a non-trivial (jr) value 
of 7 indicates the link where a dangling bond exists. 

4. Conclusion 

In conclusion, we have studied topological invariants at 
some ideal surfaces of Ge, Si described by the tight bind- 



Surface 


7 = 7T 


n e per spin 


W a 7^0 


Si, Ge, CB (111) 
Si, Ge, CB (110) 
Si, Ge (100) 
CB (100) 


71 

71,72 
71.72 

71:72 


n e = 2 
n e = 4 

n e = 4 


Wi = 1 
Wi = 1,W 4 = -1 
Wi = 1,W 2 = 1 
Wi = 1,W 2 = 1 



Table I. We have summarized material, their surface, the Berry 
phase 7 that have nontrivial value it, number of edge states n e 
per spin, and winding numbers W a that are non-zero. Here, we 
dropped a spin index because the results are spin independent: 
la = 7aT = 7aJ. an d VV Q = W a | = W a |. CB denotes the 
decoupled covalent-bond model defined in text. 



ing model with Slater and Koster's parameter. Following 
the previous study on the quantized Berry phase, 9 we 
have defined the surface-specified Berry phase. It should 
be emphasized that the symmetry required for quanti- 
zation of the Berry phase is the inversion and conjugate 
complex symmetry while it is the particle-hole symmetry 
and conjugate complex symmetry in the previous study. 
It will suggest a possibility for studies of the other ma- 
terials without the particle-hole symmetry. 

From the adiabatic transformation, Ge, Si, and the de- 
coupled covalent-bond model with the same surface are 
topologically equivalent from a viewpoint of the Berry 
phase. In other words, these topological invariants suc- 
cessfully deduced the simple picture of dangling bonds 
from the complicated tight binding model with a spin- 
orbit coupling. Moreover, the results show that the Berry 
phase is successful for all surfaces of Ge, Si and provide 
information about a position of dangling bonds. 

Finally, we comment on GaAs. For the (111) surface, 
there are four edge states n e = 4. The value of energy 
split of edge bands is roughly estimated as difference be- 
tween on site potential of two atoms. The site poten- 
tial makes it difficult the adiabatic transformation to 
the simple model as is the case with Si, Gc (100). The 
Berry phase cannot be quantized because GaAs has no 
inversion symmetry. However, the Berry phase defined 
by another twist angle can overcome the difficulty if a 
anti-unitary symmetry for it's quantization is founded in 
GaAs, 
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